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Abstract 

We present the results of our calculation of the effects of dynamical coupling of a charge-carrier 
to the electronic polarization and the field-induced lattice displacements at the gate-interface of an 
organic field-effect transistor (OFET). We find that these interactions reduce the effective band- 
width of the charge-carrier in the quasi-two dimensional channel of a pentacene transistor by a 
factor of two from its bulk value when the gate is a high-permittivity dielectric such as (Ta20s) 
while this reduction essentially vanishes using a polymer gate-insulator. These results demonstrate 
that carrier mass renormalization triggers the dielectric effects on the mobility reported recently 
in OFETs. 
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I. INTRODUCTION 



Organic field-effect transistors (OFETs) with high charge-carrier mobilities are essential 
components for high-speed organic electronic applications. For this reason, it is of crucial 
importance to unravel the origin of charge transport in these devices. Several experimental 
studies have shown that the difference in the charge-carrier transport observed in bulk or- 
ganic semiconductors and within the quasi-two dimensional channel near the gate-insulator 
interface in OFETs is associated with the effects of disorder and interfacial traps 

nag. 

More intriguing is a growing body of evidence demonstrating the strong influence of the 
dielectric permittivity of the gate- insulator on the charge- carrier mobility in OFETs j^, H| . 
Veres et al. p have studied the case of triarylamine polymer transistors in which measured 
mobilities well below 10~ 2 cm 2 /V s can be unambiguously attributed to hopping between 
localized states. On the other hand, Stassen et al. [f| obtained much higher values, from 
1 to 20 cm 2 /V s, in single- crystalline rubrene transistors in which conduction is more in- 
trinsic. A surprising finding in both cases is the drastic decrease of the mobilities as the 
dielectric constant of the gate-insulator is increased systematically using different dielec- 
tric materials. Such observation in devices governed by vastly different charge-transport 
mechanisms strongly suggests an effect due to the interactions of the charge carriers with 
the gate-dielectric. These interactions can lead to the renormalization of the bare mass of 
the charge carriers in the conducting channel of OFETs, providing an explanation for the 
observed dielectric effects on the mobility as reported in Refs. 

In this work, we present a theoretical calculation of the renormalized band mass of charge 
carriers in the conducting channel of OFETs by looking at two important interactions ex- 
perienced by charges at the interface of an organic semiconductor with a dielectric material. 
The first one, purely electronic in origin, is the image force due to the polarization discon- 
tinuity at the gate-dielectric interface. The second involves the Coulomb interaction of the 
charge carriers with surface polar phonons of the dielectric. To illustrate the basic con- 
cepts of our calculations and their general applicability to any organic semiconductor, we 
choose pentacene as a model system since its lattice structure, bandwidth, polarizabilities, 
and Huang-Rys factors are well-known compared with other organic semiconductors. We 
find that, in pentacene, the cumulative effect of the electronic and surface-lattice polaronic 
interactions is to reduce the effective bandwidth by a factor of two as the relative static 
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dielectric constant of the gate materials is varied from 1 to 25. Within the tight-binding 
model, such a decrease in the transfer integral can be viewed as a concomitant increase in 
the effective mass by the same factor. We will suggest that such renormalization of the band 
properties triors the effect, reported in Ref, gfl. 



II. TIME SCALES 



The non-interacting band properties of a perfect pentacene crystal along all crystallo- 

Q 

graphic directions have been calculated by Cheng et al. |2fl. For the bare transfer integral, 
J, between molecules along the direction of easy propagation in the (a, 6)-plane, they ob- 
tain J = 100 meV from which one gets h/J ~ 4 x 10~ 14 s as the characteristic time for 
in-plane Bloch-wave formation. The corresponding transfer time in the perpendicular c-axis 
direction, h/J±, is thirty times longer than in the plane. Thus, in the presence of scat- 
tering which substantially reduces the Bloch-wave lifetime, the carrier motion is essentially 
two- dimensional . 

The above considerations allow for the classification of the various interactions experi- 
enced by charge carriers in organic semiconductors. For fast interactions with characteristic 
times shorter than h/J, the charge can be assumed to be located on a single molecular 
site. In pentacene, this is the situation encountered during the interaction of the carrier 
with the electronic polarizability of the medium or in intramolecular charge-transfer as well 
as the coupling with intramolecular carbon stretching vibrations with frequencies around 
1360 cm -1 . Since fast interactions arise prior to the formation of the Bloch-wave, they have 
the effect of dressing the charge with a polarization cloud or a lattice deformation cloud. 
Slow interactions, on the other hand, have characteristic times much longer than h/J. They 
act directly on the Bloch-wave or the localized state. Such is the case for interaction of 
the charge carrier with low-energy intermolecular thermal phonons and librations which, in 
many cases, can be considered as static with respect to the two-dimensional band motion. 
These interactions scatter the Bloch-wave or localize the electronic states when the disorder 
they introduce is large enough. The interaction of the charge with the surface polar phonons 
of the gate insulator occurs in the intermediate time scale regime. An interesting discussion 
of time scales can also be found in the first chapter of the book by Silinsh and Capek jsj]. 

Because they dress the charge with a polarization cloud or lattice deformation, fast pro- 
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cesses lead to a renormalization of the bare transfer integrals J and J± and consequently 
increase the effective mass along all crystal directions. The case involving electron-phonon 
interactions has been discussed by several authors, including Appel ^] and Davydov |lf| . 
The purely electronic effects were treated by three of us in an earlier work ^j] in which 
we calculated the renormalization effect due to the electronic polarizability in the bulk of 
the organic semiconductor. In this work, these calculations will be extended to the situa- 
tion encountered by carriers at the gate-dielectric interface in OFETs. We shall treat both 
electronic and lattice effects. The Frohlich surface polaron at the oxide surface was already 
studied by Kirova and Bussac [12j for an isotropic organic crystal. The entire problem will 
be revisited here for the two-dimensional layer of the anisotropic crystal. The slow processes 
involving low-energy phonons, librations and other quasi-static or static sources of scattering 
and localization will be discussed and treated elsewhere. 

We have to emphasize that, because they are faster than 10 -14 s, the polarization pro- 
cesses studied here involve only the high-frequency dielectric response of the materials con- 
stituting the interface and not the usual low frequency (static) permittivity. 

For the sake of clarity and to highlight the important aspects of this work, we summarize 
the major results of our calculations in the next section. The details of these calculations 
will be given in separate appendices. 



III. RESULTS 



In rubrene- or acene-based organic field-effect transistors, the interface with the oxide or 
polymer gate-insulator involves the highly conducting (a, 6)-plane of the organic semicon- 
ductor. Our calculation of the variation of the effective transfer integral, Jjy, for conduction 
in this plane is shown in Fig^ as a function of the static dielectric constant of the gate- 
insulator. In the results of FigHJ the charge-carrier is assumed to be located on the first 
monolayer close to the gate-interface. In the case where the charge is on the second mono- 
layer, we find that there is basically no effect of the dielectric on the transfer integral which 
then takes the bulk value. These results are obtained from the combined effects of four 
different interactions which have been treated separately according to their time scales from 
the fastest to the slowest as discussed below. 

1. Electronic polarization ( J 1 } . The dynamical renormalization of the carrier motion due 
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to electronic polarization is different in the bulk and at the interface. At an interface, 
an image field is generated which is attractive when the high-frequency dielectric 
constant, eoo, of the gate- insulator is greater than that of the organic semiconductor 
as is the case for oxides, and is repulsive otherwise as in the case of polymers or air-gap 
insulators. This image field is typically of the same order of magnitude as the applied 
gate fields to which it is added or subtracted. The magnitude of the image potential 
under usual experimental conditions is displayed in FigOas a function of the distance 
z from the interface. At large distances, the classical expression for the image field 
holds and the image potential is written as, 



where €00,1 and €00,2 are the high-frequency dielectric constants of the semiconductor 
and the gate-insulator, respectively. However, corrections to this expression associated 
with lattice effects show up close to the interface. 

The increase in the carrier effective mass due to the electronic polarization cloud is 
slightly different in the bulk or when it crosses the interface. The details of these 
calculations are given in Appendix [XJ They lead to a renormalized intermolecular 
charge-transfer integral J 1 < J. 

2. Electronic displacement (J 11 ) ■ The strong dependence of the mobility and the effective 
mass with the dielectric permittivity of the gate-insulator seen in experiments, as well 
as large corrections to the bulk electronic polarization energy near the interface shown 
in FieM suggest that the first two monolayers next to the dielectric interface dominate 
the charge transport particularly in the presence of a significant gate field |2| . At even 
higher fields, one may expect that not only is the charge localized on the first monolayer 
but its electronic wavefunction is squeezed towards the part of the molecule closer to 
the insulator. This displacement of the charge distribution on the molecule is also a 
fast process, controlled by the transfer integral tn ~ 1 eV within the molecule which 
is more than an order of magnitude larger than J 1 . This fast process decreases further 
the transfer integral to J . The recent semi-empirical quantum-chemistry calculation 
performed by Sancho-Garcia et al. [l^ suggests that this effect is completely negligible 
in pentacene where the charge-carrier distribution remains perfectly centered even at 




(1) 
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very high fields (100 MV/cm). In this case Jjj = Jj. The final result of our calculation 
presented in Fig. ^ was established in the frame of this hypothesis that we consider 
as the most reliable at the moment. Nevertheless, in Appendix iBl another approach is 
presented which shows that much larger effects would be expected if the charge-carrier 
distribution on the molecule is allowed to be displaced by values of a few angstroems 
at high gate fields of 10 MV/cm. 

3. Intramolecular vibrations (J m ). Intramolecular vibrations close to 1360 cm -1 in pen- 
tacene are strongly coupled to the carrier because they change the 7r-alternation typical 
of conjugated molecules. Because these atomic motions are faster than the electronic 
polaron motion defined by the renormalization transfer integral J 11 , they also con- 
tribute to a further reduction of J 11 by a constant factor of 0.75, independent of the 
distance from the interface as well as the applied field (J 7// ~ 0.75J 77 ). Appendix O 
reviews how this factor is calculated. 

4. Frohlich polaron at the oxide surface (J IV )- Oxides are polar materials. Thus, the 
infrared-active phonon modes which modulate the metal-oxide bonds are strongly 
coupled to charge carriers sitting at their surface. In aluminum oxide, for instance, 

n 

the most active mode of this kind is situated at 46 meV [14J. This value is of the 
same order of magnitude as the effective, in-plane, transfer integral renormalized upon 
corrections due to the above-mentioned interactions. 

The construction of a lattice deformation cloud in the oxide is the object of Appendix 
IdI The polarization interaction energy of the charge with this cloud causes further 
attraction of the charge to the surface and to the subsequent increase of the effective 
mass. The calculation is performed in the intermediate coupling regime because here 
the coupling parameter which controls the process, a e s(z), defined in Appendix [Dj is 
of the order of unity in the first monolayer, and of the order of 0.1 in the second one. 

The total binding energy of the carrier in the presence of a gate-insulator includes both 
the electronic and surface polaron effects arising from the electronic image force potential 
and the lattice deformation potential at the interface associated with the above interactions. 
Within a tight-binding model, these effects are incorporated into a renormalized effective 
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transfer integral Jjv given by, 

j'\ fj"\ i j ,v 



Table HTT1 provides a summary of all these factors in the order of increasing time-scale char- 
acterizing each interaction. 

IV. CONCLUSION 

The experimental results have clearly shown that the mobilities obtained in organic field- 
effect transistors are much larger in devices built with an air-gap or a polymer gate-insulator 
□ □ than those using high-permittivity oxide gate dielectrics. Our theoretical results 
discussed above provide some insight into the origin of these effects. 

Given that the bandwidth is four times the transfer integral as suggested in the quantum 
chemistry calculations of Ref. and starting from a value of 390 meV []|, the effective 
bandwidth becomes 231 meV in bulk pentacene It is reduced to 155 meV in an OFET 
with an aluminum oxide gate insulator (AI2O3), to 146 meV close to a Ta20s interface and 
144 meV close to a TiO"2 interface, while the bulk value of 231 meV is recovered close to a 
parylene gate insulator. 

For a given disorder potential in the organic semiconductor, localization effects roughly 
scale with the reciprocal bandwidth. Consequently, important bandwidth reductions en- 
hance all localization effects and consequently decrease the mobility. It is important to note 
that in the present work the existence of a lattice polaron in the bulk of pentacene was con- 
sidered as unlikely. Intramolecular vibrations are too fast to produce a polaron in a perfect 
pentacene crystal and their effect has been studied in Appendix Intermolecular phonons 
and librations with energies fkv of the order of 10 meV are not coupled enough to the car- 
rier to produce a well-defined lattice polaron, with a renormalized transfer integral of the 
order of 50 meV originating from the initial calculation of Cheng et al. |7j. However, when 
the transfer integral experiences a further reduction close to the gate interface as shown in 
Fig. [TJ then the intrinsic lattice polaron can be excited, as suggested in a recent work jlfi| . 
A calculation is in progress to clarify this point. Moreover, an interface is rarely perfect 
from a structural point of view, and the fact that, due to significant internal field present 
at the interface, the carrier is constrained to probe more closely all the interface disorder, 
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enhancing the localization effects. Even "perfect" interfaces can be "intrinsic" sources of 
localization. In general, the electric dipole lattice induced by a charge carrier in the organic 
semiconductor is incommensurate with the dipole lattice induced in the gate material. In 
a perfect structure, the incommensurability of the electronic potential can open gaps in 
the semiconductor density of states. When the gate dielectric is not a single- crystal but a 
disordered structure, the same effect creates an electronic disorder which can be one of the 
sources of localization. 

The effective bandwidth enters the mobility law but does not define it entirely. In disor- 
dered polymer channels, Veres et al. [^| have established a link between the effective width 
of the Gaussian distribution of electronic states and the mobility of the carrier which jumps 
from state to state according to the Gaussian Disorder Model [Ijj • In more ordered systems 
such as single- crystalline rubrene OFETs, a theoretical work is in progress to establish a 
quantitative link between the transfer integral and the mobility. It will elucidate the role of 
thermal, low-energy phonons which are, on the one side, sources of localization and on the 
other the key to the adiabatic diffusion of the carriers in the channel. 
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APPENDIX A: ELECTRONIC POLARON AT THE INTERFACE 

The electronic polarization of a molecular crystal, i.e. the formation of induced dipoles 
d on the neutral molecules of the crystal in the field of an extra-charge, occurs on a time 
scale of r ~ 10~ 15 — 10~ 16 s (h/r ~ 2 eV). This is much faster than the carrier lifetime on 
a molecule which is typically T/~/i/J~4x 10 -14 s with J ~ 100 meV, allowing one to 
treat the carrier as stationary on a molecular site during the polarization process. 
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1. Bulk 



With the above approximation, we have calculated the polarization energy of a point- 
charge in bulk pentacene crystal to be E p ~ —1.5 eV ^l] using the method of self-consistent 
polarization field to the dipolar approximation 8]. The pentacene crystal structure and pa- 
rameters are given in Tab HI The molecules involved in the discrete calculation are assumed 
to be inside a spherical cluster of radius R. The polarization energy varies linearly with 
the inverse radius of the cluster. The polarizabilities of the molecules along the (L, M, N)- 
directions were adjusted so that the polarization energy at large R matches the continuum 
limit (l — Crff 1 ) e 2 /47reo-R, where e e fr is the effective dielectric constant in an anisotropic struc- 
ture jlSj. The convergence is fast and the calculated polarization energies for anthracene, 
E r = -L26 eV, and pentacene, E r * -1.5 eV, are in good agreement with recently reported 
values |19j . 

As the charge carrier moves from one molecule to another, the polarization energy 
changes. In fact, when the electronic polarization energy and the induced-dipole distri- 
bution vary from site n to site n + h due to different local crystal environment, the carrier 
will encounter a resistance to its motion. This leads to a renormalization of the bare transfer 
integral J. We introduced the parameter S |nj as a way of describing quantitatively the 
polarizability effect on the transfer integral, i.e. J 1 = Jexp (—So). The quantity S is closely 
related to the relative amplitudes and directions of the induced dipoles and is given through 
the relations [ll|. 

1 3 

SoW = ^EE (") - X M (" + h )f ( A1 ) 

i 1=1 

Xil = (A2) 

where an are the components of the polarizability, and e, of the order of one electron-volt, 
is the energy difference between the ground state and the first-excited state of the neutral 
molecule The sites considered here, n and n + h, are the next-nearest neighbor along 
the crystal direction d 2 = — + |b. 
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2. Interface 



We now extend the above calculation to the case of charge carrier near the interface 
of pentacene with the dielectric insulator. In this case, the polarization cloud extends 
partly into the semiconductor and partly into the dielectric. For simplicity, we modelled 
the dielectric insulator as a cubic lattice. The Clausius-Mosotti relation is then applied to 
obtain the electronic polarizability of the molecules of the dielectric assuming a given high 
frequency permittivity, 

eoo,2 - 1 1 



£00,2 + 2 3eo 



Na e (A3) 



N = ^N A (A4) 

Here M and p are the molar mass and the density of the dielectric and Na is the Avogadro 
number. The summation in Eq JAll is then extended to include the induced dipoles in the 
semiconductor and the dielectric. For reasons of symmetry, the molecules involved in the 
discrete calculations are taken inside a cylindrical cluster on both sides of the interface. 
Beyond the cluster approximation, a continuum contribution is added to the polarization 
energy using standard electrostatics 20]. The polarization energy as a function of the 
distance from the interface is shown in FigJ2]for both the discrete and the continuous limits. 



APPENDIX B: SQUEEZING OF THE CHARGE DISTRIBUTION AT THE IN- 
TERFACE 

The self-trapping interfacial polarization electric field (see Appendix EJ and the applied 
gate-field localize the excess charge on the pentacene molecule by pulling and squeezing the 
wavefunction towards or away from the interface depending on the relative magnitudes of 
the dielectric constants of pentacene and the gate- insulator. To get the average position and 
the spatial extent of the wavefunction of the excess charge, we used a model in which the 
pentacene molecule is considered as a pair of conjugated chains of eleven sites, each separated 
by a distance / that we treat in the tight-binding approximation. The basis wavefunction 
on each site, ip n , with n — 1, .., 11 satisfies the following set of equations, 

E ip n = -t// (V>n+i + VVi-i) - nlqF?p n - V (nl) i/; n (Bl) 
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with the boundary conditions ijj = ^12 = imposed at the two ends of the pentacene 
molecule. Here V is the image force potential, F is the applied gate-field, and tu is the in- 
tramolecular transfer integral. Upon solving this system of equations, we obtain the ground- 
state wavefunction as, 

11 



|*o> = 5>», (B2) 



i=l 



The spatial dependence of the average position of the excess charge on the chain will be 
given by, 

11 

i) = Y^nm 2 (B3) 



11 

71=1 

Since the pentacene molecule is nearly perpendicular to the interface, the presence of 



interfacial fields localizes the charge at the extremity of the molecule. This has the effect 
of changing the thickness of the channel and, as a consequence, the bandwidth for charge 
propagation. Indeed, the transfer integral in acene crystals increases with the number of 
aromatic rings from naphtalene to pentacene 0, 1^- This can be related to the spatial 
extent of the excess charge on the acene molecule. The charge extent a n is defined as the 
width containing 80% of the charge density. We have chosen as in Appendix El the direction 
of higher transfer integral which determines the bandwidth. Therefore, the variation of 
the charge extension under the effect of the image force and the gate-field, will induce a 
corresponding change in the transfer integral of the pentacene molecules as shown in Fig. 
El For the calculation of the excess charge extension in the above model, the intramolecular 
transfer integral t u was set to leV. The values of the gate and image force fields are those 
found in Appendix El However, a semi-empirical quantum chemistry calculation published 
recently has shown that, probably due to exchange interactions between the carrier and 
the 7r-electrons in the pentacene molecule, the extra-charge distribution is not pulled and 
squeezed in the direction of the interface field by the extent predicted by the above model. 
Thus in the final result of Fig. ^ this effect has not been included. Figure El shows what 
would the final result be in the case where the squeezing calculated above in this appendix 
would indeed be effective. 
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APPENDIX C: BAND NARROWING DUE TO INTRAMOLECULAR VIBRA- 
TIONS 



We now refer to the strong coupling of the extra charge to the intramolecular vibrations 
at 1360 cm -1 in pentacene. The coupling constant has been determined in acenes both 



221. T 



experimentally and theoretically in Ref. 
phonons has been presented by several authors 
for dispersionless intramolecular vibrations. 



re calculation of the band narrowing due to 
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101 ] . Here, we present a simplified version 



The electron-phonon Hamiltonian of interest is written as, 

H = ^ -J Ha n+h a n + X] hUoKbn ~ \/E B /hhjQ (&+ + b n ) d+t 



(CI) 



n,h 



where n represents the molecular sites; a+, a n , b+, and b n are the electron and phonon 
operators respectively, g = y^E B /fku is the usual coupling parameter. Since the phonon 
frequency huio is larger than the transfer integral J 11 , we look for a variational solution of 
the Hamiltonian. The trial wavefunction is of the form, 



where 



l^n) = ^U n \n) <g> \Xn) 



| X „)=exp (X*b n -X n b+)\0) 



(C2) 



(C3) 



describes the intramolecular vibrations of the molecule n on which the charge is located. 
The variationnal parameter X n is determined by minimizing the energy 



e = mm 



(C4) 



Then 



cxp 



n,h 



(- 




2 + 


X n+ h 


2 \ 


\ 2 ; 



U n+h U n 



-y/E B /?ku (X n + X*) \u n \ 2 + ^o\X n 

n 

We find X n = ^E B /huo and the corresponding transfer integral, 



(C5) 



J 1U = J u exp(-E B /huJc 



(C6) 
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When thermal phonons are taken into account the thermal average value of the energy yields 
a temperature correction factor in the transfer integral given by, 



j "'= jn ^{-tM^f)) (C7) 

In the present case, the temperature effect is completely negligible and the electron-phonon 
binding energy is just — Eb per molecular site. 

Taking the values in bulk pentacene from Ref. with Eb = 45 meV, we obtain 

J 111 = 0.75J 77 . Close to the interface the band narrowing due to intramolecular vibrations 
is rather insensitive to the interfacial field. This is related to the fact that the phonon 
frequency and coupling constants are approximately the same in the acene series |22j |. 

APPENDIX D: THE FROHLICH SURFACE POLARON 

When a charge carrier is generated by the field effect at the interface between a molecular 
semiconductor and a gate insulator,it interacts with the surface phonons of the dielectric. 
This effect has been studied in Ref. |l2| for an isotropic 3D molecular crystal in the adiabatic 
limit. Here, we consider this interaction in the case of a pentacene crystal where the carrier 
motion is essentially two-dimensional and for moderate electron-phonon coupling. Indeed, 
the residence time in a monolayer r ~ h/J±, with J± ~ 5 meV is much larger than the time 
to polarize the dielectric given by 27r/u s , where hw s = 46 meV in AI2O3 jl4j . However, the 
phonon frequency Ttw s is of the same order of magnitude as the effective in-plane transfer 
integral J 111 ~ 60 meV. The frequency uj s is given by the following formula: 

^ = 1k + 4), (di) 



where ujl and ujt are the frequencies of the bulk longitudinal and transverse phonons. [23j 
The values of uj s for the different dielectrics are reported in Tab. [H] 

The electron-phonon interaction involving a charge in a particular monolayer of the crystal 
at a distance z > from the interface is given by |8[ , 



q » n 



where \ip n \ is the charge density at site na = (n x a, n y a). Here b q and b_ q are the annihilation 
and creation operators of the surface phonons in the gate material, u s their frequency, and 
S the surface area of the interface. 
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If e x (w) is the dielectric susceptibility of the molecular crystal, €00 2 and e ,2 the high and 
low frequency limits of the dielectric permittivity, the coupling constant 1/e* is given by, 



1 



ei — £oo, 2 



e* ei (ex + 600,2) ei (ei + e , 2 ) 
and the total Hamiltonian becomes, 



(D3) 



e-ph 



(D4) 



where p is the momentum of the charge carrier. Following Ref. 



24j . we introduce the total 



momentum of the system which is a constant of motion of the total Hamiltonian, 



p = h ^K h * + p 



(D5) 



We can transform the total Hamiltonian H to H' through the unitary transformation S, so 
that H' no longer contains the charge coordinates, 



H' = S^HS 



with, 



We obtain thus, 



S = exp 



(D6) 



(D7) 



h' = y, ^Kb q +J2 V « (*) + b - q ) 



P/h-J2bp q q 



J in a 



+ 



1IIJ1 



J J11 a 



where 



V q {z) = — 



-qz 



(D8) 



(D9) 



v^V Se* 

With the phonon frequency hw s comparable to the effective transfer integral J 111 , the adi- 
abatic approximation is not applicable. However, the dimensionless parameter a e g which 
describes the strength of the electron-phonon coupling decreases with the distance z to the 
interface as 

e 2 exp (— 2irz/a) 



"off 



87re e*a ^hw s J lu 



ae 



-2-irz/a 



(D10) 
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In our case, a e g is of the order of 1 for a charge carrier located in the first monolayer. Then, 
we use a variational method to describe the interaction of the dressed charge carrier with 
the dielectric phono. Introducing a second unitary —.nation, 



i> = exp^6 J +/,-6,/ I «j 



where /„ will be chosen to minimize the energy 



(Dll) 



E = -a'J lu + J2 (V q f q + V q *f* q ) + J 111 £ | /g |Va 2 



h 



twj. + J 1 " UV-2q~a 2 

fr 



We then find, 



fq 



V* 



lq huj s + J in 

where 77 satisfies the implicit equation 

77P 



q 2 a 2 - ^Y~a 2 (1 - 77) 



(D12) 



(D13) 



(D14) 



hw s + J IH [q 2 a 2 - 2^-a 2 (1 - rj)] 

The carrier binding energy is obtained as Ef, = —al\ (z) ftw s and the effective mass is 
m*/m = J IH /J IV = 1 + 20/2(2). As long as J IH P 2 a 2 /h 2 is small (< Huj s ), we may 
obtain E (P 2 ) to first order in an expansion in powers of Oj^^fif-^j ■ On doing so, one 



readily gets 



where 



and 



2 ? 



E = —ali (z) fkj s + 



P*a 



J 



in 



h 2 [l + 2al 2 (z) 



h 0) 



h(z) 



"y/J" 1 /*" dy 



1 + y- 



exp I ~7V ^ > 9 ^ 



(l + y 2 f eXP [-^\j^ y]giy) 



Here the term given by 



g(y) 



sinh 



,, , huj s e 



, huj s t 



(D15) 



(D16) 



(D17) 



(D18) 
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accounts for the finite extension of the charge distribution on the molecule of length t. 



[1] V. Podzorov, E. Menard, A. Borissov, V. Kiryukhin, J. A. Rogers, and M. E. Gershenson, 

Phy. Rev. Lett. 93, 086602 (2004). 
[2] F. Dinelli, M. Murgia, P. Levy, M. Cavalleni, F. Biscarini, and D. M. de Leeuw, Phys. Rev. 

Lett. 92, 116802 (2004). 

[3] J. H. Kang, D. A. da Silva Filho, J. L. Bredas, and X.-Y. Zhu, Appl. Phys. Lett. 86, 152115 
(2005). 

[4] M. Daraktchiev, A. von Miihlenen, F. Niiesch, M. Schaer, M. Brinkmann, M. N. Bussac, and 

L. Zuppiroli, New Journal of Physics 7, 133 (2005). 
[5] J. Veres, S. D. Ogier, S. W. Leeming, D. C. Cupertino, and S. M. Khaffaf, Adv. Func. Mater. 

13, 199 (2003). 

[6] A. F. Stassen, R. W. I. de Boer, N. N. Iosad, and A. F. Morpurgo, Appl. Phys. Lett. 85, 3899 
(2004). 

[7] Y. C. Cheng, R. J. Silbey, D. A. da Silva Filho, J. P. Calbert, J. Cornil, J. L. Bredas, J. Chem. 

Phys. 118, 3764 (2003). 
[8] E. A. Silinsh and V. Capek, Organic Molecular Crystals (AIP Press, New York, 1994). 
[9] J. Appel, Solid Stat. Phys. 25, 193 (1968). 
[10] A. S. Davydov, Theorie du Solide (Edition Mir, 1980). 

[11] M. N. Bussac, J. D. Picon, and L. Zuppiroli, Europhys. Lett. 66, 392 (2004). 
[12] N. Kirova, and M. N. Bussac, Phys. Rev. B 68, 235312 (2003). 

[13] J. C. Sancho-Garcia, G. Horowitz, J. L. Bredas, and J. Cornil, J. Chem. Phys. 119, 12563 

(2003) . 

[14] M. Schubert, T. E. Tiwald, and C. M. Herzinger, Phys. Rev. B 61, 8187 (2000). 
[15] D. A. da Silva Filho, E. G. Kim, and J. L. Bredas, Adv. Mater. 17, 1072 (2005). 
[16] K. Hannewald, and P. A. Bobbert, Phys. Rev. B 69, 075212 (2004); K. Hannewald, V. M. 
Stojanovic, J. M. Schellekens, P. A. Bobbert, G. Kresse, and J. Hafner, ibid. 69, 075211 

(2004) . 

[17] H. Bassler, Phys. Stat. Sol. B 175, 11 (1993). 

[18] P. J. Bounds, and R. W. Munn, Chem. Phys. 44, 103 (1979). 

16 



[19] E. V. Tsiper, and Z. G. Soos, Phys. Rev. B 68, 085301 (2003). 
[20] J. D. Jackson, Classical Electrodynamics (Wiley & Sons, New York, 1999), p. 154. 
[21] J. L. Bredas, J. P. Calbert, D. A. da Silva Filho, and J. Cornil, Proc. Natl. Acad. Sci. 99, 
5804 (2002). 

[22] V. Coropceanu, M. Malagoli, D. A. da Silva Filho, N. E. Gruhn, T. G. Bill, and J. L. Bredas, 

Phys. Rev. Lett. 89, 275503 (2002). 
[23] J. Sak, Phys. Rev. D 6, 3981 (1972). 

[24] T. D. Lee, F. E. Low, and D. Pines, Phys. Rev. 90, 297 (1953). 



17 



TABLE AND FIGURE CAPTIONS 



Figure 1: The largest transfer integral for charge propagation in pentacene close to 
the interface with the gate insulator as a function of the relative static permittivity 
of the dielectric. The charge carrier is located on the first monolayer close to the 
interface. The high frequency dielectric constant €00 2 for the gate insulators is: 2.65 
for parylene C; 2.37 for Si0 2 ; 3.1 for A1 2 3 ; 4.12 for Ta 2 5 ; and 6.5 for Ti0 2 . 

Figure 2: The image potential E p (z) at the interface is given in the continuous 
approximation and in our lattice model. At large distances from the interface the 
electronic polarization energy in the bulk is recovered. 

Figure 3: These are the same results as in Fig. ^but calculated with another set of 
data as explained in Appendix El The largest transfer integral for charge propagation 
in pentacene close to the interface with the gate insulator is plotted as a function of 
the relative static permittivity of the dielectric. The charge-carrier is assumed to be 
located on the first monolayer. The inset depicts the case where the charge- carrier is 
on the second monolayer. These curves cumulate all the effects calculated in this work. 
In contrast to the results of Fig. E]we have assumed here that the carrier distribution is 
pulled away or squeezed towards the interface according to the calculation in App endix 
iBl This effect is enhanced with respect to the semi-empirical model of Ref . |lj| . Here 
the effective mass depends even on the gate field. At present, the corresponding values 
calculated in Fig. dare considered more reliable. 

Table 1: Crystal constants of pentacene [R. B. Campbell et al., Acta Cryst. 14, 705 
(1961)]. 

Table 2: Values of the surface phonon frequencies used in the calculations of Appendix 
iDl together with the coupling constant a (see text). 

Table 3: Reduction factors of the transfer integral through the steps discussed in 
the text. Cases of bulk pentacene and interfaces with vacuum and different dielectric 
oxides are shown along with the time scales characterizing each process. Here the 
second effect which is related to charge displacement on the molecule was excluded 
(see text). 
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Crystal constants (A) 


a 





c 




f .y 


O.UO 


10. Ul 


V_/l_ysLrdl CUIi.oLaJ.lLo ^ULg^ 




P 


7 




101.9 


112.6 


85.8 


Dielectric constant 0. 


L 


M 


AT 




5.336 


3.211 


2.413 



a E.V. Tsiper and Z.G. Soos, Phys. Rev. B 68, 085301 (2003). 

TABLE I: 





uj s (cm 1 ) 


a 


Si0 2 a 


480 


2.84 


A1203 6 


386 


5.24 


Ta 2 5 c 


390 


6.17 


Ti0 2 d 


280 


4.84 



a J. Humli'cck, A. Roseler, Thin Solid Films 234, 332 (1993) 
b M. Schubert, ct al., Phys. Rev. B 61, 8187 (2000) 
C E. Franke, et al., J. Appl. Phys. 88, 5166 (2000) 
d R. Sikora, J. Phys. Chcm. Solids 66, 1069 (2005) 

TABLE II: 





Molecular 
polarization 


Intramolecular 
charge vibration 


Surface phonons 


final result 


Typical time scales 


~ 10- 15 s 


~2x 10- 14 s 


~8x 10- 14 s 




Reduction factors 


J 1 /J 


J 111 /J 1 


J IV /.J UI 


J IV /J 


Pentacene bulk 


0.79 


0.75 


1 


0.593 


Pentacene /vacuum 


0.76 


0.75 


1 


0.57 


Pentacene /SiO 2 


0.77 


0.75 


0.815 


0.47 


Pentacene/Al2C>3 


0.77 


0.75 


0.69 


0.398 


Pentaccnc/Ta20s 


0.77 


0.75 


0.65 


0.375 


Pcntaccnc/Ti02 


0.77 


0.75 


0.64 


0.37 



TABLE III: 
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Bulk value 
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FIG. 3: 
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